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The LiH + molecule is prototypical of the indirect dissociative recombination (DR) process, in 
which a colliding electron destroys the molecule through Rydberg capture pathways. This Letter 
develops the first quantitative test of the Siegert state multichannel quantum defect theory descrip- 
tion of indirect DR for a diatomic molecular ion. The R-matrix approach is adopted to calculate 
ab-initio quantum defects, functions of the internuclear distance that characterize both Rydberg 
states and the zero-energy collisions of electrons with LiH + ions. The calculated DR rate coefficient 
agrees accurately with recent experimental data We identify the doorways to fast indirect DR 
as complex resonance manifolds, which couple closed channels having both high and low princi- 
pal quantum numbers. This sheds new light on the competition between direct and indirect DR 
pathways, and suggests the reason why previous theory underestimated the DR rate by an order of 
magnitude. 

PACS numbers: 



The LiH molecule was one of the first participants in 
early universe chemistry, and it played a key role in the 
cooling of primordial gases. Therefore, the relative abun- 
dance of LiH and its formation process have drawn sig- 
nificant attention in models of the universe. Stancil et 
al. @ noticed that direct radiative association of neutral 
atoms, Li + H — > LiH + v occurs at the very slow rate 
coefficient of ~ 1CP 20 cm s . 

On the other hand, radiative association with ionic 
hydrogen Li + H + — > LiH + + v is predicted to occur 
d, [1] at the much higher rate of ~ 10~ 15 cm 3 s -1 . The 
resulting LiH + and LiH abundances are controlled by 
photoionization, collisions, and by dissociative recombi- 
nation (DR) with free electrons. A recent DR experi- 
ment for LiH+ Q has measured the DR rate coefficient 
for collisions at T = 139K (12 meV) to be (6±2) x 10~ 7 
cm 3 s _1 . From numerous ab-initio calculations of the LiH 
and LiH + potential surfaces [f| H, 0, H, 0, EH there is no 
Born-Oppenheimer neutral state that crosses the ionic 
ground state potential curve, anywhere near the ionic 
minimum. Thus the high measured DR rate for such an 
indirect or non-crossing process is challenging to recon- 
cile with existing theoretical results. The rate coefficient 
estimated theoretically in Q is 2.6 x 10 -8 cm 3 s _1 . An- 
other theoretical study by Florescu et al. applied 
multichannel quantum defect methods, with the relevant 
nonadiabatic coupling elements obtained from a gener- 
alized Hellmann-Feynman theorem 11| to calculate that 
the 139K DR rate coefficient should equal 3.6 x 10~ 8 
cm 3 s _1 . Since both of these theoretical studies underes- 
timate the DR rate for this simple diatomic by more than 
an order of magnitude, it shows that the proper physi- 
cal description of the underlying mechanism for indirect 
DR processes continues to challenge our theoretical un- 
derstanding. 

Quantum defect theory in connection with frame trans- 
formation into a basis of Siegert vibrational states 12, HJ 



has shown promise in describing DR for a model diatomic 
[3] and for the triatomic H^" molecule systems for 
which indirect Rydberg state pathways dominate. But to 
date there has been no rigorous test of the Siegert-state- 
based MQDT formulation, for an experimentally-studied 
system for which the relevant quantum defect matrices 
have been determined directly in an ab initio scattering- 
type calculation. Accordingly, the main goal of the cur- 
rent paper is to analyze dissociative collisions between a 
low-energy electron and the LiH + ion, as a fundamen- 
tal prototype system that provides a stringent test of 
this combination of theoretical elements: R-matrix the- 
ory, MQDT, and frame-transformation theory based on 
Siegert pseudostates 1(3, 17, 3, 19j |. 



The body-fixed adiabatic eigenquantum defects n^{R) 
are related to the energy differences (in a.u.) between 
the potential curves of the ion U + (R) and the neutral 
Rydberg states U^AR) by Mulliken's equation 



U+(R) 



2[<(i?)]' 
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where the effective quantum number v^(R) = n — n^{R) 
and A denotes the projection of the Rydberg electron an- 
gular momentum I onto the axis of the diatomic molecule. 
Of course I is not a good quantum number, so 7 is 
an eigen-index distinguishing different Rydberg series of 
LiH. Fundamentally, body-frame quantum defects are 
represented by a matrix /i^, (R) with different partial- 
waves coupled by off-diagonal elements, after which the 
(J>~(R) are obtained as its eigenvalues. Ab-initio poten- 
tial energy surfaces by themselves provide no direct in- 
formation about off-diagonal couplings, or equivalently, 
about the eigenvectors C// 7 (i?) , and many authors tend to 
neglect them (T3] or estimate them using a two-channel 
Landau-Zener curve-crossing formula [201 ] or sometimes 



by fitting them in a diabatic representation 211 ]. 
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But the present study exploits the familiar MQDT 
theorem [18| that smoothly connects quantum defects 
(multiplied by tt) at energies just below the ionization 
threshold to i?-dependent short-range scattering phase 
shifts (or multichannel scattering or reaction matrices) 
just above the threshold: 



tion and boundary conditions: 



it 4,{R) = J2 U h (R) 6$(R) U Vl {R) 



(2) 



Here 6^ are the low- or zero-energy eigenphases for e - 
+ LiH + collision and the eigenvector matrix is U(R) 
which transforms the short-range if-matrix into diago- 
nal form. We used the diatomic UK R-matrix package 
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FIG. 1: Upper panels: The calculated quantum-defect matrix 
elements are shown versus the internuclear distance - panel 
a) singlets, panel b) triplets. Lower panels: The solid curves 
are the eigenvalues of the matrices in a) and b). The black 
dots are eigenquantum defects extracted from n=4 states in 
extensive bound state CI calculations [ljj. 

[22T | to calculate the short-range if-matrix by matching to 
coulomb functions at an R-matrix boundary of ro — 25ag. 
The target was described by an augVTZ STO basis set 
[23I ]. From the calculated 1 S + and 3 E + quantum de- 
fects shown in Fig. [TJ the s-wave quantum defect is only 
weakly perturbed from its Li + limit of 0.399, over a wide 
range of internuclear distances. Thus coupling to nuclear 
motion will be controlled by higher partial waves, namely 
p- and d-waves. The 1 H and 3 II quantum defects have 
been calculated to be an order of magnitude smaller, and 
they have negligible impact on the final DR results. 

The internuclear distance R is a good body frame 
"quantum number" when all electrons are confined 
within the box specified by ro and the Born-Oppenheimer 
approximation is strictly valid. The vibrational frame 
transformation connects R with the laboratory-frame 
quantization of this degree of freedom expressed by vi- 
brational wave functions <pj V (R). Siegert pseudostates 
[lil [u| provide a unified description of the bound vi- 
brational states and the vibrational continuum. In the 
present case they solve the vibrational Schrodinger equa- 



d 2 i(i + 1) , 2 

3 V 
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4>jv(o) = o 
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R 2 
d_ 
dR 



4> jv (R)=0,(3) 



- ikj V (j) JV (R) 



= 0.(4) 



Ro 



In the above equations j is a rotational quantum num- 
ber of the ion, M stands for its reduced mass, while 
Ro denotes a nuclear radius beyond which we approxi- 
mate the interaction potential in ([3]) to be constant and 
4>j V (R) = exp(ikj V R) for R > Rq. Fig. [2] shows an 
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FIG. 2: Distribution of the j = Siegert momentum eigen- 
values from Eqs. (|3l4 |l in the complex plane. The circled states 
are the ones used in the calculations of this paper. 

example of the j — Siegert state momentum eigen- 
value distribution. Because the nuclei are confined within 
Rq = 10ao we only obtain 4 bound states, in contrast 
with the expected total of 7 bound states found in @, [l(| ■ 
However for the lower vibrational levels that fit inside, 
agreement in the level spacing is achieved within 2 cm . 

Because the orthogonality relation between two differ- 
ent Siegert pseudostates is slightly modified a sur- 
face term [14| is added to the standard frame transfor- 
mation integral [24[ yielding 



Si,vAh?) = j Q dR^R)^^)^^) 



+ 1- 



(5) 



The underline in this equation denotes that fi A is a ma- 
trix with indices ^{R). Moreover, the rotational in- 
dices (j,j') of this body-frame 5-matrix do not give 
rotational transition probabilities - they only serve as 
a reminder that the vibrational functions exhibit a j- 
dependence through the centrifugal term in Eq.([3]). To 
reiterate, the rotational frame transformation transforms 
a set of body-frame quantum numbers (I, A, J) into a set 
of laboratory-frame quantum numbers (l,j,J). The to- 
tal angular momentum J is defined via J = I + j. The 
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LiH + ion is treated in Hund's case (b), with spin-orbit 
coupling neglected. Definite total parity 77 = (— is 
also enforced, whereby the short-range laboratory-frame 
scattering matrix is: 



qJt) 

ljv,V j'v' 



(0) 



The real, 



Ujjl is taken from 



orthogonal rotational transformation matrix 
ill ] and will not be repeated here. 
As is familiar in MQDT applications, the "short- 
range" or "unphysical" scattering matrix S_ in Eq.© is 
diagonal in the J and rj quantum numbers. It represents 
an amplitude for electron-ion scattering from an initial 
channel defined by (I'j'v') into a final channel defined 
by (Ijv), but some of these channels in S_ are typically 
closed energetically. The physically-relevant S'-matrix is 
of course defined only in the open-channel space and is 
obtained by the "elimination of closed channels" formula 
Iz|: 



gphys _ go 



S' c 



2i0(E) 



e - 



S c 



(7) 



where the superscripts o and c respectively denote open 
and closed sub-blocks of the unphysical S'-matrix ^ and 
(3(E) is a diagonal matrix of effective Rydberg quantum 
numbers with respect to the closed-channels thresholds: 



ft- 



irdi 



V 2 ( E i - E ) 



(8) 



The total energy of the electron + ion system is E, and 
Ei is the ionization threshold for channel i = (vjl). Here 



as in Refs. [14 , |15| | the high ionization thresholds are de- 
scribed by a Siegert pseudo-continuum state with finite 
widths, so Ei and (3(E) are complex. This fact alone de- 
stroys the unitarity of S phys making it sub-unitary. The 
lost flux is associated with a trapped Rydberg electron in 
a closed channel that represents a high-lying vibrational 
state that is dissociative, with outgoing-wave character, 
and has a complex vibrational energy and corresponding 
finite lifetime. The departure from unitarity was identi- 
fied by 14, 1|| as the dissociation probability following 



electron impact in incident channel i': 



i-E^ ys (£)4f ys (£) 



(9) 



with the incident electron collision energy ey = E—Eji v i . 
This cross-section depends on the initial channel i' — 
(v'j'l'); the collision preserves the conserved quantum 
numbers {Jrj). The experimentally observable cross- 
section for dissociation following electron impact is then 



Tji v i ( 



(10) 



We further average over a Boltzmann distribution of ini- 
tial ro-vibrational states of the ion at the temperature 
(T = 300K) appropriate to the experiment [H. The cor- 
responding recombination rate a(E c ) is 



t{E c ) = y^2E c a(E c ) 



(11) 



This calculated DR rate exhibits an infinite number of 
resonances near each closed-channel ionization threshold, 
associated with autoionizing and predissociating states of 
LiH. To compare with the storage ring experiments [l[, 
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FIG. 3: DR rate: The solid curve is our calculated, 
anisotropically-averaged rate for AE«= 0.1 meV and AE± 
= 12 meV. The broken curve shows the calculation of fiol ] . 
The dotted curve is our truncated result obtained by neglect- 
ing the off-diagonal couplings in the quantum-defect matrix 
(see the fig. [TJ. The crosses denote the experimental data 
[261 ]. with a few characteristic error bars shown. 

we must convolve over an anisotropic finite spread in the 
electron energy; the spread is different for the parallel 
(A£j| =0.1 meV) and the much broader AE±^ = 12 meV 
perpendicular components of the relative velocity vector. 
The convolution over parallel and perpendicular energy 



rjjl' 



distributions has been performed as was outlined in [15 1 
and elaborated in detail in [251 ] . Figure [3] summarizes 
our results along with previous experimental and theo- 
retical results. This figure also demonstrates the results 
of a numerical test conducted to interpret the discrep- 
ancy between our theoretical results and those of [101 ] . 
Specifically, we have performed one set of calculations 
that neglect the off-dia gon al Z-mixing to mimic the cal- 
culations performed by [101 ] . i.e. using only the diagonal 
eigenvalue form of quantum-defect matrix (shown 

in fig. [T]c,d). Introduction of this approximation lowers 
the DR rate by an order of magnitude, and this arti- 
ficially restricted calculation agrees with the results of 
[1(1 ]. Thus, the presence of i?-dependent Z-mixing is cru- 
cial for this system. Moreover, the rate is predominantly 
controlled by (i-wave collisions, whereas Ref.[l3| assumed 
that it was dominated by p- waves. 

We now identify the qualitative mechanism responsi- 
ble for this high indirect DR rate. Figure 3] summarizes 
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the probabilities of various DR pathways at an energy 
near the first vibrational threshold. For clarity we have 
neglected rotations in this qualitative analysis, because 
we found very little effect of the rotational frame trans- 
formation on the present results at this energy resolu- 
tion. In FigH]no thermal averaging has been applied, and 
the LiH + ion is initially in the vibrational ground state. 
The DR probabilities shown are the quantity inside the 
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FIG. 4: Contributions to DR probabilities are shown as func- 
tions of energy near the first vibrational threshold: Tht thick 
broken curve represents the average DR probability. The 
full spectrum of Rydberg resonances converging to the first 
threshold is denoted by the thin line. Contributions from 
the second and higher thresholds are shown as the thick solid 
curve. The thin broken curve is the probability of a vibra- 
tional excitation process vo — ► which is only energetically 
allowed above the v = 1 threshold. 

square brackets in Eq.©. The thick broken curve de- 
notes the average DR probability across the threshold. 
The DR probability drops from about 5% of the inci- 
dent flux below the threshold down to only 0.6% above 
the threshold. The DR probability below the threshold 
is built up as the cumulative effect of the dense forest 
of Rydberg resonances attached to the first vibrationally 
excited state, seen as the thin full curve. The thick full 
curve shows the contributions to DR probability when 
the incident electron is captured into Rydberg states as- 
sociated with higher vibrational thresholds. FigfJ] also 
explains the reason why the DR flux drops sharply above 
the v = 1 threshold. The thin broken curve shows our 
calculated probability of vibrational excitation — ► 1. 
As can be noticed from the amplitudes for the probabil- 
ities of both processes below and above the threshold, 
most of the DR flux just below the threshold turns dis- 
continuously into vibrational excitation flux once that 
channel becomes open. These results indicate that the 
DR process is controlled by a doorway, namely capture 
of the incident electron into a Rydberg state attached to 
the first vibrational threshold. However, if there is no 
higher- v/lower-n perturbing level overlapping the total 
energy of the collision complex, the electron will tend to 
autoionize before DR can take place. Throughout the en- 
ergy range of a multichannel complex resonance, though, 



the initial capture can efficiently pump more energy into 
vibration at the first electron recollision, and Figf4]shows 
that this resonant perturbed Rydberg complex increases 
the DR rate by about another factor of 3. This complex 
resonance mechanism for indirect DR is believed to ap- 
ply to many other systems that are not controlled by the 
usual simple capture mechanism into a dissociative state. 
A hint of its importance in H-j" DR is evident across the 
energy range 110-160 cm -1 in Fig. 10 of Ref . [27I] . 
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